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CLAIMS 

[1] A signal separation method that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 
5 a procedure that transforms the observed signal values observed by 

said sensors into frequency-domain signal values, 

a procedure that uses said frequency-domain signal values to calculate 
at each frequency the relative values of the observed values between said 
sensors (including mapping these relative values), 
10 a procedure that clusters said relative values into N clusters, 

a procedure that calculates a representative value for each of said 
clusters, 

a procedure that uses said representative values to generate a mask for 
the purpose of extracting, from said frequency-domain signal values, mixed 
15 signal values comprising the signals emitted from V (2<V<M) signal sources, 

a procedure that uses said mask to extract said mixed signal values 
from said frequency-domain signal values, and 

a procedure that separates and extracts the values of V signals from 
said mixed signal values. 
20 [2] A signal separation method according to Claim 1 , 

wherein said mask is a function that takes a high level value for said 
relative values that are within a prescribed range that includes V said 
representative values, and takes a low level value for said representative 
values that are not inside said prescribed range, 
25 and wherein the procedure that uses said mask to extract said mixed 

signal values from said frequency-domain signal values is a procedure in 
which said frequency-domain signal values are multiplied by said mask. 
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[3] A signal separation method according to Claim 1, 

wherein said mask is a function that takes a low level value for said 
relative values that are within a prescribed range that includes V said 
representative values, and takes a high level value for said representative 
5 values that are not inside said prescribed range, 

and wherein the procedure that uses said mask to extract said mixed 
signal values from said frequency-domain signal values is a procedure in 
which the values obtained by multiplying said frequency-domain signal 
values by said mask are subtracted from said frequency-domain signal values. 

10 [4] A signal separation method according to Claim 2, 

wherein said mask is a function that the transitions from said high 
level value to said low level value that accompany changes of said relative 
value occur in a continuous fashion. 
[5] A signal separation method according to Claim 1, 

1 5 wherein the procedure that uses said representative values to generate 

a mask for the purpose of extracting, from said frequency-domain signal 
values, mixed signal values comprising the signals emitted from V (2<V<M) 
signal sources is a procedure whereby said mask is generated by using the 
directional characteristics of a null beamformer (NBF). 

20 [6] A signal separation method according to Claim 1 , wherein the 
procedure that uses said representative values to generate a mask for the 
purpose of extracting, from said frequency-domain signal values, mixed 
signal values comprising the signals emitted from V (2<V<M) signal sources 
includes: 

25 a procedure that generates an (N-V+l)x(N— V+l) delay matrix 

H NB F(f) in which the element at (j,i) is equal to exp(j27ifiji), where 
Tji=(dj/v)cos0i, v is the velocity of the signals, dj is the distance between 
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sensor 1 and sensor j (j=l ,...,N-V+1), 0i is any one of the estimated directions 
of the signal sources corresponding to the V said representative values, 0i 
(i=2,...,N-V+l) are the estimated directions of the signal sources 
corresponding to the other said representative values of the V said 
representative values, and f is a frequency variable, 

a procedure that calculates the inverse matrix W(f)=H NB F" 1 (f) of delay 
matrix H NB F(f) as a NBF matrix W(f), 

a procedure that generates a directional characteristics function 
FORMULA 54 

N-V+l 



10 F(f,0)= £ W lk (f)exp(j27cfd k cos0/v) 



k=l 



where 0 is a signal arrival direction variable, and the first row element of said 
NBF matrix W(f) is W lk (f), 

and a procedure that uses said directional characteristics function 
F(f,0) to generate said mask. 
15 [7] A signal separation method according to Claim 1 , wherein the 
procedure that uses said representative values to generate a mask for the 
purpose of extracting, from said frequency-domain signal values, mixed 
signal values comprising the signals emitted from V (2<V<M) signal sources 
includes: 

20 a procedure that generates a function consisting of a single-peak 

function convolved with a binary mask, which is a function that takes a high 
level value for said relative values that are within a prescribed range including 
V said representative values and takes a low level value for said 
representative values that are not inside said prescribed range and where 

25 changes of the relative value are accompanied by discontinuous transitions 
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from said high level value to said low level value, 

and a procedure that generates said mask as a function in which said 
relative values are substituted into said function consisting of a single-peak 
function convolved with a binary mask. 
5 [8] A signal separation method according to Claim 1, wherein the 
procedure that uses said representative values to generate a mask for the 
purpose of extracting, from said frequency-domain signal values, mixed 
signal values comprising the signals emitted from V (2<=V<=M) signal 
sources is 

10 a procedure that generates said mask as a single-peak function 

obtained by mapping the differences between a first odd function that takes a 
value of zero when said relative value is the lower limit value a min in a 
prescribed range including V said representative values and a second odd 
function that takes a value of zero when said representative value is the upper 

15 limit value a max in said prescribed range. 

[9] A signal separation method according to Claim 2 or Claim 3, wherein 

said mask is a function that transitions from said high level value to 
said low level value occur in a discontinuous fashion 

[10] A signal separation method that separates and extracts signals under 
20 conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 

a procedure that transforms the observed signal values observed by 
said sensors into frequency-domain signal values, 

a procedure that uses said frequency-domain signal values to calculate 
25 at each frequency the relative values of the observed values between said 
sensors (including mapping these relative values), 

a procedure that clusters said relative values into N clusters, 
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a procedure that calculates a representative value for each of said 
clusters, 

a procedure that generates a mask function that takes a high level 
value for said relative values that are within a prescribed range that includes 
5 one of the said representative values, and takes a low level value for said 
representative values that are not inside said prescribed range, wherein the 
transitions from said high level value to said low level value that accompany 
changes of said relative value occur in a continuous fashion, 

and a procedure that multiplies said frequency-domain signal values 
10 by said mask to extract the signal emitted from one signal source. 

[11] A signal separation method that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 

a procedure that transforms the observed signal values Xi(t),...,x M (t) 
1 5 observed by said sensors into frequency-domain signal values 
X i (f,m), . . . ,X M (f,m), 

a procedure that clusters first vectors X(f,m)=[X 1 (f,m),...,X M (f,m)] 
comprising said frequency-domain signal values Xi(f,m),...,X M (f>m) into N 
clusters Q(f) (i=l,...,N) at each frequency f, 
20 a procedure that calculates second vectors aj(f) to represent each said 

cluster Ci(f), 

a procedure that extracts V (1<V<M) third vectors a p (f) (p=l,...,V) 
from said second vectors aj(f), 

a procedure that generates a mask M(f,m) represented by the formula 
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FORMULA 55 

M(f,m) = j 1 maX VOec k D(X(f,m),a p (f)) < ™\ {{)eGt D(X(f,m),a q (f)) 
[O otherwise 

where G k is the set of said third vectors a p (f), G k ° is the complementary set of 
G k , and D(a,{3) is the Mahanalobis square distance between the vectors a and 
5 (3, 

and a procedure that extracts the signal values emitted from V of said 
signal sources by calculating the product of said mask M(f,m) and said first 
vectors X(f,m). 

[12] A signal separation method that separates and extracts signals under 
10 conditions where N (N>2) signals are mixed together and observed with M 
sensors, wherein 

a procedure that transforms the observed signal values Xi(t),...,x M (t) 
observed by said sensors into frequency-domain signal values 
X j (f,m), . . . ,X M (f ,m), 

15 a procedure that clusters first vectors X(f ? m)=[Xi(f,m) 3 ...,X M (f,m)] 

comprising said frequency-domain signal values Xi(f,m),...,X M (f ? m) into N 
clusters Q(f) (i=l,...,N) at each frequency f, 

a procedure that calculates second vectors a;(f) to represent each said 
cluster Cj(f), 

20 a procedure that extracts V (1<V<M) third vectors a p (f) (p=l,...,V) 

from said second vectors ai(f), 

and a procedure that judges whether or not said first vectors X(f,m) 
satisfy the relationship 
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max 



Vf)eG k 



D(X(f,m),a p (f)) < min aq(f)eG , D(X(f,m),a q (f)) 



where G k is the set of said third vectors a p (f), G k c is the complementary set of 
Gk, and D(a,(3) is the Mahanalobis square distance between the vectors a and 
5 (3, and, if so, extracts said first vectors X(f,m) as the signal values emitted 
from V of the said signal sources. 

[13] A signal separation method according to Claim 1 1 or Claim 12, 
wherein said clustering procedure is performed after performing the 
calculation 
10 FORMULA 57 



[14] A signal separation method according to Claim 13, wherein said 
1 5 clustering procedure is performed after performing the calculation 




(| Xjftrn)!* 0) 
(| XjftnOH 0) 



and 
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X(f,m)<- 




(where the notation || X(f,m) || denotes the norm of X(f,m)). 



After said formula 
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fXff.mVsienfX. 

X(f,m)<- 



X(f,m)/sign(X j (f,m» (| X^f.m)^ 0) 
X(f,m) (1^(^111)1=0) 



[15] A signal separation method that separates and extracts signals 
under conditions where N (N>2) signals are mixed together and observed 
with M sensors, comprising 
5 a procedure that transforms the observed signal values Xi(t),...,x M (t) 

observed by said sensors into frequency-domain signal values 
X i (f,m), . . . ,X M (f>m), 

a procedure that clusters first vectors X(f,m)=[Xi(f,m),...,X M (f,m)] 
comprising said frequency-domain signal values Xi(f,m),...,X M (f,m) into N 
10 clusters Ci(f) (i=l,...,N) at each frequency f, 

a procedure that calculates second vectors ai(f) to represent each said 
cluster Q(f), 

a procedure that calculates an N-row x M-column separation matrix 
W(f,m) that is the Moore-Penrose pseudo-inverse matrix of an M-row x N- 
15 column matrix in which 0 or more of the N said second vectors ai(f) are 
substituted with zero vectors, 

and a procedure that calculates a separated signal vector 
Y(f,m)=[Y 1 (f,m) 5 ...,Y N (f,m)] by performing the calculation 
Y(f,m)=W(f,m)X(f,m). 
20 [16] A signal separation method according to Claim 15, wherein: 

the procedure that calculates said separation matrix W(f,m) is a 
procedure that selects min(M,N) said second vectors aj(f), generates a matrix 
A'(f,ni) whose columns are the selected min(M,N) said second vectors aj(f) 
and max(N-M,0) zero vectors, and calculates said separation matrix W(f,m) 
25 as the Moore-Penrose pseudo-inverse matrix of said matrix A'(f>m). 
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[17] A signal separation method according to Claim 15, wherein: 

the procedure used to calculate said separation matrix W(f,m) when 
N>M is a procedure that selects M said second vectors aj(f) in each discrete 
time interval m, generates a matrix A'(f,m) whose columns are the selected M 
said second vectors aj(f) and N-M zero vectors, and calculates said (time- 
dependent) separation matrix W(f,m) as the Moore-Penrose pseudo-inverse 
matrix of said matrix A'(f,m), 

and the procedure used to calculate said separation matrix W(f,m) 
when N<M is a procedure that calculates the Moore-Penrose pseudo-inverse 
matrix of a matrix comprising N said second vectors in each said cluster Q(f) 
to yield said (time-invariant) separation matrix W(f,m). 
[18] A signal separation method according to Claim 15, wherein: 

said clustering procedure is performed after performing the calculation 



[19] A signal separation method according to Claim 18, wherein: 

said clustering procedure is performed after performing the calculation 
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X(f,m)<- 
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(where the notation || X(f,m) || denotes the norm of X(f,m)). 
in addition to said formula 



X(f,m) 



'X(f, m) / signer (f, m)) (| Xj(f, m) |* 0) 
X(f,m) (|X j (f,m)|=0) 



[20] A signal separation method according to Claim 16, wherein 
5 said procedure that selects min(M,N) said second vectors ai(f) is a 

procedure that initializes fourth vectors e with said first vectors X(f,m), and 
then repeats a process min(M,N) times in which it selects said second vectors 
a q ( u) (f) that maximize the absolute value of the dot product of a q ( U )(f)/||a q(U )(f)|| 
and said fourth vectors, sets up a matrix Q = [a q( i)(f) ? ... ? a q ( k )(f)] representing the 
10 subspace subtended by all said second vectors a q(u )(u=l v ..,k) selected so far, 
performs the calculation P=Q(Q H Q) _1 Q H , and updates the fourth vectors e 
with e=X(f,m)-P-X(f,m). 

[21] A signal separation device that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
15 sensors, comprising: 

a memory unit that stores the observed signal values observed by said 
sensors; 

and a processor which is connected to said memory unit and performs 
processing whereby it; 
20 transforms said observed signal values into frequency-domain signal 

values, 

uses said frequency-domain signal values to calculate at each 
frequency the relative values of the observed values between said sensors 
(including mapping these relative values), 
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clusters said relative values into N clusters, 

calculates a representative value for each of said clusters, 

uses said representative values to generate a mask for the purpose of 

extracting, from said frequency-domain signal values, mixed signal values 
5 comprising the signals emitted from V (2<V<M) signal sources, 

uses said mask to extract said mixed signal values from said 

frequency-domain signal values, and 

separates and extracts the values of V signals from said mixed signal 

values. 

10 [22] A signal separation device that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 

a memory unit that stores the observed signal values observed by said 
sensors; 

15 and a processor which is connected to said memory unit and performs 

processing whereby it; 

transforms said observed signal values into frequency-domain signal 

values, 

uses said frequency-domain signal values to calculate at each 
20 frequency the relative values of the observed values between said sensors 
(including mapping these relative values), 

clusters said relative values into N clusters, 
calculates a representative value for each of said clusters, 
generates a mask, which is a function that takes a high level value for 
25 said relative values that are within a prescribed range that includes one said 
representative value, and takes a low level value for said representative values 
that are not inside said prescribed range, and where the transitions from said 
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high level value to said low level value that accompany changes of said 
relative value occur in a continuous fashion, 

and extracts the values of a signal emitted from one signal source by 
multiplying said frequency-domain values by said mask. 
5 [23] A signal separation device that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: . 

a memory unit that stores the observed signal values xi(t),...,x M (t) 
observed by said sensors; 
10 and a processor which is connected to said memory unit and performs 

processing whereby it; 

transforms said observed signal values Xi(t),...,x M (t) into frequency- 
domain signal values Xi(f,m),...,X M (f 9 m), 

clusters first vectors X(f,m)=[Xi(f,m),...,X M (f 5 m)] comprising said 
15 frequency-domain signal values Xi(f,m),...,X M (f,m) into N clusters Q(f) 
(i=l,...,N) at each frequency f, 

calculates second vectors aj(f) to represent each said cluster Q(f), and 
extracts V (1<V<M) third vectors a p (f) (p=l,...,V) from said second vectors 



ai(f), 



20 



generates a mask M(f,m) represented by the formula 
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1 maX a p (f)eG k 



otherwise 




25 



where G k is the set of said third vectors a p (f), G k ° is the complementary set of 
G k , and D(a,(3) is the Mahanalobis square distance between the vectors a and 

p 
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and extracts the signal values emitted from V of the said signal 
sources by calculating the product of said mask M(f,m) and said first vectors 
X(f,m). 

[24] A signal separation device that separates and extracts signals under 
5 conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 

a memory unit that stores the observed signal values Xi(t),...,x M (t) 
observed by said sensors; 

and a processor which is connected to said memory unit and performs 
1 0 processing whereby it; 

transforms said observed signal values Xi(t),...,x M (t) into frequency- 
domain signal values Xi(f,m),... 5 X M (f,m), 

clusters first vectors X(f 5 m)=[Xi(f,m) v .. ? X M (f,m)] comprising said 
frequency-domain signal values X 1 (f,m),...,X M (f,m) into N clusters Ci(f) 
15 (i= 1 , . . . ,N) at each frequency f, 

calculates second vectors aj(f) to represent each said cluster Ci(f), 

extracts V (1<V<M) third vectors a p (f) (p=l,...,V) from said second 
vectors aj(f), 

judges whether or not said first vectors satisfy the relationship 
20 FORMULA 62 

max a p( f)eG k D ( x ( f ' m )^ a p( f ))<min aq(f)eG , D(X(f,m),a q (f)) 

where G k is the set of said third vectors a p (f), G k c is the complementary set of 
Gk, and D(a,f3) is the Mahanalobis square distance between the vectors a and 
P, and, if so, extracts said first vectors X(f,m) as the signal values emitted 
25 from V of the said signal sources. 
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[25] A signal separation device that separates and extracts signals under 
conditions where N (N>2) signals are mixed together and observed with M 
sensors, comprising: 

a memory unit that stores the observed signal values Xi(t),...,x M (t) 
5 observed by said sensors; 

and a processor which is connected to said memory unit and performs 
processing whereby it; 

transforms said observed signal values Xi(t),...,x M (t) into frequency- 
domain signal values Xi(f,m),...,X M (f,m), 
10 clusters first vectors X(f 5 m)=[Xi(f,m),...,X M (f,m)] T comprising said 

frequency-domain signal values Xi(f,m),...,X M (f,m) into N clusters C\(f) 
(i=l,...,N) at each frequency £, 

calculates second vectors ai(f) to represent each said cluster Ci(f), 

calculates an N-row x M-column separation matrix W(f,m) that is the 
1 5 Moore-Penrose pseudo-inverse matrix of an M-row x N-column matrix in 
which 0 or more of the N said second vectors aj(f) are substituted with zero 
vectors, 

and calculates a separated signal vector Y(f,m)=[Y 1 (f,m),... ? Y N (f ? m)] T 
by performing the calculation Y(f,m)=W(f,m)X(f,m). 
20 [26] A signal separation program that causes a computer to perform: 

a procedure that transforms observed signal values, which are 
mixtures of N (N>2) signals observed with M sensors, into frequency-domain 
values, 

a procedure that uses said frequency-domain signal values to calculate 
25 at each frequency the relative values of the observed values between said 
sensors (including mapping these relative values), 

a procedure that clusters said relative values into N clusters, 
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a procedure that calculates a representative value for each of said 
clusters, 

a procedure that uses said representative values to generate a mask for 
the purpose of extracting, from said frequency-domain signal values, mixed 
5 signal values comprising the signals emitted from V (2<V<M) signal sources, 

a procedure that uses said mask to extract said mixed signal values 
from said frequency-domain signal values, and 

a procedure that separates and extracts the values of V signals from 
said mixed signal values. 
10 [27] K signal separation program that causes a computer to perform: 

a procedure that transforms observed signal values, which are 
mixtures of N (N>2) signals observed with M sensors, into frequency-domain 
values, 

a procedure that uses said frequency-domain signal values to calculate 
15 at each frequency the relative values of the observed values between said 

sensors (including mapping these relative values), 

a procedure that clusters said relative values into N clusters, 

a procedure that calculates a representative value for each of said 

clusters, 

20 a procedure that generates a mask, which is a function that takes a 

high level value for said relative values that are within a prescribed range that 
includes one of said representative values, and takes a low level value for said 
representative values that are not inside said prescribed range, wherein the 
transitions from said high level value to said low level value that accompany 

25 changes of said relative value occur in a continuous fashion, 

and a procedure that extracts the signal values emitted from one signal 
source by multiplying said frequency-domain signal values by said mask. 
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[28] A signal separation program that causes a computer to perform: 
a procedure that transforms observed signal values xi(t),...,x M (t), 

which are mixtures of N (N>2) signals observed by M sensors, into 

frequency-domain signal values X 1 (f,m),...,X M (f,m), 

a procedure that clusters first vectors X(f,m)=[Xi(f,m) v ..,X M (f,m)] 

comprising said frequency-domain signal values Xi(f,m),...,X M (f,m) into N 

clusters Ci(f) (i=l v ..,N) at each frequency f, 

a procedure that calculates second vectors ai(f) to represent each said 

cluster Q(f), 

a procedure that extracts V (1<V<M) third vectors a p (f) (p=l,...,V) 
from said second vectors ai(f), 

a procedure that generates a mask M(f>m) represented by the formula 
FORMULA 63 



where G k is the set of said third vectors a p (f), G k ° is the complementary set of 
Gk, and D(a,(3) is the Mahanalobis square distance between the vectors a and 



and a procedure that extracts the signal values emitted from V of said 
signal sources by calculating the product of said mask M(f,m) and said first 
vectors X(f,m). 

[29] A signal separation program that causes a computer to perform: 
a procedure that transforms observed signal values Xi(t),...,x M (t), 

which are mixtures of N (N>2) signals observed by M sensors, into 

frequency-domain signal values Xi(f,m),...,X M (f,m), 

a procedure that clusters first vectors X(f,m)=[Xi(f,m),...,X M (f,m)] 




0 otherwise 
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comprising said frequency-domain signal values Xi(f,m),...,X M (f,m) into N 
clusters C\(f) (i=l v ..,N) at each frequency f, 

a procedure that calculates second vectors a\(f) to represent each said 
cluster Ci(f), 

a procedure that extracts V (1<V<M) third vectors a p (f) (p=l,...,V) 
from said second vectors aj(f), 

and a procedure that judges, whether or not said first vectors X(f,m) 
satisfy the relationship 
FORMULA 64 

max a p( f)eG k D ( x ( f ' m X a p( f ))<min aq(f)eG , D(X(f,m),a q (f)) 

where Gk is the set of said third vectors a p (f), G k c is the complementary set of 
G k , and D(a,p) is the Mahanalobis square distance between the vectors a and 
p, and, if so, extracts said first vectors X(f,m) as the signal values emitted 
from V of the said signal sources. 

[30] A signal separation program that causes a computer to perform: 
a procedure that transforms observed signal values X!(t),...,x M (t), 

which are mixtures of N (N>2) signals observed by M sensors, into 

frequency-domain signal values Xi(f,m),...,X M (f,m), 

a procedure that clusters first vectors X(f,m)=[X 1 (f,m),...,X M (f,m)] 

comprising said frequency-domain signal values Xi(f,m),...,X M (f,ni) into N 

clusters Q(f) (i=l,...,N) at each frequency f, 

a procedure that calculates second vectors aj(f) to represent each said 

cluster Ci(f), 

a procedure that calculates an N-row x M-column separation matrix 
W(f,m) that is the Moore-Penrose pseudo-inverse matrix of an M-row x N- 
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column matrix in which 0 or more of the N said second vectors aj(f) are 
substituted with zero vectors, 

and a procedure that calculates a separated signal vector 
Y(f,m)=[Yi(f 5 m) v ..,Y N (f,m)] T by performing the calculation 
Y(f,m)=W(f,m)X(f,m). 

[31] A computer-readable recording medium that stores a signal separation 
program according to any one of Claims 26 through 30. 



